Finite-temperature properties of quasi-2D Bose-Einstein condensates 
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Using the finite-temperature path integral Monte Carlo method, we investigate dilute, trapped 
Bose gases in a quasi-two dimensional geometry. The quantum particles have short-range, s-wave 
interactions described by a hard-sphere potential whose core radius equals its corresponding scatter- 
ing length. The effect of both the temperature and the interparticle interaction on the equilibrium 
properties such as the total energy, the density profile, and the superfluid fraction is discussed. We 
compare our accurate results with both the semi-classical approximation and the exact results of 
an ideal Bose gas. Our results show that for repulsive interactions, (i) the minimum value of the 
aspect ratio, where the system starts to behave quasi-two dimensionally, increases as the two-body 
interaction strength increases, (ii) the superfluid fraction for a quasi-2D Bose gas is distinctly dif- 
ferent from that for both a quasi-lD Bose gas and a true 3D system, i.e., the superfluid fraction for 
a quasi-2D Bose gas decreases faster than that for a quasi-lD system and a true 3D system with 
increasing temperature, and shows a stronger dependence on the interaction strength, (iii) the su- 
perfluid fraction for a quasi-2D Bose gas lies well below the values calculated from the semi-classical 
approximation, and (iv) the Kosterlitz-Thouless transition temperature decreases as the strength of 
the interaction increases. 
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The study of the effects of reduced dimensionality 
has attracted considerable interest in condensed mat- 
ter physics and statistical physics, for example, quantum 
films of superfluid '^He on surfaces, superfluid ''He clus- 
ters, and superfluid *He confined in restricted geometries. 
Recently, the experimental realization of Bose-Einstein 
condensation (BEC) in dilute, trapped, and supercooled 
atomic vapors 1] and the additional ability to produce a 
low dimensional atomic gas trapped in an optical lattice 
have stimulated the investigation in a low dimensional 
system both theoretically 13, p, Ij , | a , Jyl and experimen- 
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I21I l22l|. Dimensionally reduced systems have very differ- 
ent properties from their three dimensional (3D) counter- 
parts 13 due to the enhanced importance of phase 
fluctuations. For example, in a spatially homogeneous 
infinite system a one-dimensional (ID) Bose gas does 
not exhibit BEC, even at zero temperature; and in two 
dimensions (2D) BEC exists only at zero temperature. 
Nonetheless, a dilute two-dimensional Bose gas under- 
goes a superfluid phase transition at a finite critical tem- 
perature. Below the Kosterlitz-Thouless transition tem- 
perature 25], the gas is superfluid, and the superfluid 
phase is characterized by the presence of a quasiconden- 
sate, i.e., a condensate with only local phase coherence, 
since long wavelength phase fluctuations destroy long- 
range order. 

Especially^ intriguing are two-dimensional Bose 
gases m IH El El S m m m m, and re- 
cently, experiments have entered regimes of BEC in 
2D||ESE1E1E1E1- In a spatially homogeneous 
two-dimensional system Safonov et al. observed 
the first experimental evidence for a quasicondensate 
in 2D atomic hydrogen gas adsorbed on a superfluid 
''He surface. Furthermore, BEC in trapped gases is 
qualitatively different from BEC in a homogeneous 



potential. Lower dimensional trapped Bose gases have 
been realized with strong quantum confinement in 
one or more directions |3(l | (i) by gradually reducing 
atoms from a highly anisotropic trap to decrease the 
interaction energy [3, Il0j |: (ii) by gradually increasing 
the trap anisotropy while keeping the number of atoms 
fixed [H E3; and (iii) by using the periodic potential 
of a one-dimensional optical lattice 0, 0| . When the 
system is in a harmonic trap, the effect of the trap 
becomes more dramatic with lowering of dimensionality 
of the system; the external trapping potential limits 
the size of the atomic gas, and the density of states is 
modified. As a result, Bagnato and Kleppner |33 showed 
that an ideal 2D gas in a harmonic potential does exhibit 
a BEC phase at a finite temperature. Petrov et al. |3^ 
showed that well below Tc the equilibrium state is a true 
condensate, whereas at intermediate temperatures a 
quasicondensate forms when local two-body interactions 
are included. 

Theoretically, the ground-state properties of a trapped 
Bose gas in two dimensions have been studied recently 
using the Gross-Pitaevskii mean-field theory |3 Issf . 
the leadin g qu antum corrections to the Gross-Pitaevskii 
equation j27| . and a variational model based on a 
Gaussian-parabolic trial wave function [23|- For the 
finite-temperature properties of a dilute Bose gas con- 
fined in a harmonic trap, the semiclassical approxima- 
tion js^, a microscopic mean field theory that includes 
both density and phase fluctuations of the Bose gas 
the Hartre-Fock-Bogoliubov formalism [26l l28l E3 j and 
the scaling structure within diagrammatic perturbation 
theory [sj have been used. 

However, many properties have yet to be investigated 
both experimentally and theoretically. Needless to say, a 
more thorough theoretical understanding is clearly desir- 
able. At finite temperatures among the main problems 



2 



to be answered are the density profile as a function of 
temperature and the strength of the interparticle inter- 
action, under what conditions the confinement gives a 
system with dimensionality 2, the superfluidity, and the 
effect of interactions on the critical temperature. In ad- 
dition, as the role of correlations and of quantum fluctua- 
tions is enhanced in quasi-2D geometries, an appropriate 
theoretical description requires the more accurate many- 
body approach beyond the mean-field approximation for 
a more systematic investigation in a wider range of tem- 
perature and values of the scattering length. 

In this paper we use a finite-temperature path-integral 
Monte Carlo (PIMC) method |^ to investigate the 
crossover from 3D to 2D for N particles, where the inter- 
particle interaction is a purely repulsive hard-sphere po- 
tential of radius as, the s-wave scattering length. PIMC 
allows one to calculate accurate quantum mechanical ex- 
pectation values of many-body system, the only input 
being the many-body potential. In particular, we demon- 
strate the infiuence of temperature T and interparticle 
interaction Os on the equilibrium properties of ultracold 
atomic gases to calculate their energetics and structural 
properties in a highly anisotropic trap. Furthermore, we 
compare our accurate PIMC results with both the semi- 
classical approximation and the exact values for an ideal 
Bose gas, and we provide detailed predictions for com- 
parison with future experiments. 

We consider the following Hamiltonian for a system of 
N hard spheres 



N 

H = Ho + J2v{n,) (1) 

i>j 
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where Hq is the Hamiltonian for trapped ideal Bose gases 
and v{r) is a two-body, symmetric hard-sphere potential 
defined by 

v(r) - / +°° < 
^ > " \ (r > as). 

In order to investigate the crossover from 3D to 2D, we 
consider Bose gases under a cylindrically harmonic con- 
finement [ujx — tOy = LiJp and Uz = AciJp), where A is the 
aspect ratio. Here we shall always be concerned with the 
properties of trapped Bose gases at finite temperatures 
with the number of particles N , the scattering length a^, 
of the two-body interaction potential, and two charac- 
teristic lengths az = yjh/muiz and Op = ^hjmujp^ de- 
scribing the oscillation lengths in the transverse and the 
longitudinal directions, respectively. We gradually vary 
the axial trap frequency to increase the trap anisotropy 
and to enter the quasi-two-dimensional regime, where the 
particles obey 2D statistics but interact in the same way 
as in a three-dimensional system. As a result, the gas 
is confined by an extreme, pancake-shaped potential and 



FIG. 1: The calculated density profiles n(p), normalized such 
that n{p)pdp = 1 and p = y/x^ + y^, at T = 0.4 Tc for 
two different aspect ratios, i.e., A = 10 and A — 100. We used 
three different scattering lengths as = (solid lines), as — 
aitb (dotted lines), and as = lOaRt (dashed lines). The inset 
shows the expectation value of p in unit of ap as a function of 
A for as — (triangles), as = ajib (diamonds), and Os = lOttRi, 
(circles). In all figures, when statistical errors cannot be seen 
on the scale of the figure, the error bars are smaller than the 
symbol sizes. The expectation value of p clearly depends on 
both the scattering length as and the aspect ratio A. The 
dependence on as becomes increasingly large as A increases. 

has quasi-two-dimensional properties. Along the tightly 
confined axial direction the motion is frozen out, so that 
the condensate is in the harmonic oscillator ground state 
of an ideal gas, and the condensate width equals the har- 
monic oscillator length. The tight confining direction is 
characterized by the frequency uj^ and two weak confin- 
ing directions by ujp. 

In this study, we typically use N—27 hard spheres be- 
cause the permutation sampling in PIMC is efficient if 
iV = L'^, where L is an integer |4j|, and the next largest 
number that satisfies this criterion would have required 
an excessive amount of cpu time. There is a finite-size 
effect in the density profile and the total energy; how- 
ever, the superfluid fraction shows no dependence upon 
N within the statistical error of the data. Experimen- 
tally, in an optical lattice, one can control and reduce the 
number of particle in each BEC 0,E3- Consequently, it 
will be possible to compare the effects predicted here for 
small numbers of atoms with experimental data directly. 
We also use several trap aspect ratios 1 < A = LUz/tOp < 
300 and a wide temperature range, i.e., 0.1 < T/Tc < 1.4, 
wheret Tc is the transition temperature for an ideal Bose 
gas of N atoms in a trap (see, e.g., Eq. (19) of Ref. 0). 
Due to the finite-size corrections, the critical temperature 
Tc depends on both the number of atoms iV and trap- 
ping frequencies, ujp and Uz- For a repulsively interacting 
trapped Bose gas, Tc is lowered compared to the trapped 
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FIG. 2: The calculated density profiles n(p) at the aspect 
ratio A = 100 for two different temperatures, i.e., T /Tc — 0.2 
and T/Tc — 1.2 using three different scattering lengths as — 
0, Us = ajit, and as = 10a_Ri,. The inset shows the dependence 
of the expectation values of p on the temperature T/Tc- The 
arrow on the y-axis indicates the expectation value of p at T 
= 0, < p > — 0.886 Op. As the temperature increases, the 
density profiles n{p) expand along the p axis and the effect of 
the interaction strength decreases. 



ideal Bose s result of interaction, in contrast to an 

uniform Bose gas in the dilute range, where the critical 
temperature is increased above the ideal Bose gas value 
by interaction. Here, the length unit is Up, where Op = 
y^h/mujp, and energies are measured in units of hcop. 

In addition, we assume that the s-wave scattering 
length aRb of ^'^Rb is 100 times the Bohr radius, i.e., 
am = 0.00433 Op. PIMC procedure used here is based 
on methods described in our previous work |43l |. 

We first discuss the density profiles of iV = 27 hard 
spheres at T = 0.4 Tc as a function of A, 1 < A < 300. 
Figure 1 shows the calculated density profiles n{p) as 
a function of p, normalized such that n{p)pdp = 1 

and p = + for two different aspect ratios, i.e., 
A = 10 and A — 100. To study the effect of the interac- 
tion strength, we used three different scattering lengths 
Os = (solid lines), as = aj^b (dotted lines), and — 
lOaRb (dashed lines). Clearly, the total density profile 
along the p coordinate spreads out in the trap as both 
as and A increases, i.e., the density at the center of the 
trap decreases with increasing A and as as expected. In 
particular, the density profile n(p) depends strongly on 
the interaction strength as for A > 10 and changes dra- 
matically as a function of the aspect ratio A. However, 
for A = 1 (a spherically symmetric harmonic trap), the 
density profile n{p) for Og = a^b approaches that of the 
non-interacting gas (see the inset of Fig. 1). 

The inset of Fig. 1 shows the expectation value of p in 
unit of flp as a function of A for as = (triangles), as 



= aRb (diamonds), and as = lOaRb (circles). In all fig- 
ures, when statistical errors cannot be seen on the scale 
of the figure, the error bars are smaller than the sym- 
bol sizes. The expectation value of p clearly depends 
on both the scattering length and the aspect ratio A. 
The dependence on as becomes increasingly large as A 
increases, which suggests that the repulsion between the 
atoms spreads the atoms in the trap as a result of inter- 
action. 

Next, we further examine the condition to achieve a 
quasi-2D Bose gas using the density profiles n{z) along 
the axial direction as a function of z, normalized such 
that J^^n{z)dz — 1, and the expectation value of \z\. 
In the quasi-2D regime, along the z direction the gas 
at finite temperatures has the characteristics of an ideal 
non-interacting gas at T = 0. Thus, the density profile 
n(z) becomes identical to that of an ideal gas, 

n{z) = ^^exp[-{z/a,f]. (2) 

From Eq. the expectation value of \z\ is 0.564 Oz for 
an ideal Bose gas at T = 0. 

Our calculated expectation values of |z| for = anb at 
T = QATc approach the ground state value (T — 0) as the 
aspect ratio A increases. Our calculated density profile 
n(z) for as = anb at T = 0.4rc and that for an ideal gas 
(Eq. (0)) arc indistinguishable for A > 40. This suggests 
that the motion along the z coordinate is indeed frozen 
out as A approaches large values since the energy of the 
axial excitations increases with increasing A and it is not 
easy for particles to go to the excited states for large val- 
ues of A. For repulsive interactions, the minimum value 
of A, where the system starts to enter a quasi-2D regime, 
increases as the interaction strength increases. 

To see the effect of the temperature T, we also cal- 
culated the density profiles n{p) and n(z) and the ex- 
pectation values of p and \z\ as a function of the scaled 
temperature T/Tc for the aspect ratio A = 100 and A'' = 
27. 

Figure 2 shows the calculated density profiles n{p) as 
a function of p, for two different temperatures, i.e., T/Tc 
— 0.2 and T/Tc = 1-2 using three different scattering 
lengths as = 0, = aj^b, and as = lOaRb (using the 
same symbols as in Fig. 1). As the temperature increases, 
the density profiles n{p) expand along the p axis and the 
effect of the interaction decreases. 

The inset of Fig. 2 shows the dependence of the ex- 
pectation values of p on the temperature T/Tc- At very 
low temperatures, our calculated expectation values of p 
for fls = approach the ground state value, i.e., {p) = 
0.886 ftp (see the arrow on the y-axis). As the tempera- 
ture increases, the expectation values of p increase, but 
the effect of the interaction decreases. The expectation 
values of p for Og = and as — a^b are the same within 
errorbars at high temperatures. However, the expecta- 
tion values of p for as = lOaub are still larger than those 
for fls = even at high temperatures. 
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FIG. 3: The calculated PIMC total energies per particle E/N 
in units of TtWp as a function of the scaled temperature T /Tc 
for N — 27 and A = 100. We use three different interaction 
strengths as — 0, as — aRt, and as = 10a jn,. In the figure, we 
subtract Eg = hujz/2 = 50 huip, the ground energy of an ideal 
Bose gas in the tight confinement direction, from our PIMC 
total energy. The arrow on the y-axis indicates the ideal gas 
energy per particle in the p direction, hup, at T — 0. 



FIG. 4: The superfluid fraction along the axis of rotation z 
for T/Tc = 0.4 and A'^ =27 as a function of the aspect ratio 
A calculated using the PIMC method. In this calculation 
we have used three different interaction strengths as = 0, 
as = aRb, and as — lOaRt- At the aspect ratios A < 100, 
the superfiuid fraction decreases as the aspect ratio increases. 
However, the superfluid fractions do not depend on the aspect 
ratio, to within error bars, for aspect ratios A > 100. 



In contrast to the expectation value of p, the expecta- 
tion value of |z| for A = 100 is nearly constant for T < 
lATc, the whole temperature range used here, and is the 
same as that for the ground state of an ideal Bose gas. 
Thus, the motion in the axial direction is largely frozen 
out and the system behaves quasi-two dimensionally at 
those temperatures. This demonstrates clearly that for 
the axial axis the system at finite temperatures behaves 
like an ideal gas at T = 0. 

Figure 3 shows our calculated PIMC total energies per 
particle E/N in units of hujp as a function of the scaled 
temperature T/T^ for iV = 27 and A = 100. We used 
three different interaction strengths as = 0, as = aflb, and 
Os = lOaRb- For large anisotropics and small atom num- 
bers, the total energy per particle tends towards Tiujz/'^, 
the ground energy of an ideal Bose gas in the tight con- 
finement direction. As shown above, the motion along 
the z direction is indeed frozen out at A = 100 for all 
temperature range used in this calculation 0.1 < T /Tc < 
1.4, i.e., the total energy per particle in the tight confine- 
ment direction equals that for an ideal Bose gas in one 
dimension, hiOz/'^ = 50 Tiu^p. In the figure, we subtracted 
i?Q = h,ijJz/'2. = 50 hijjp from our PIMC total energy per 
particle E/N to interpret our results. 

Our calculated PIMC energies per particle for a = 
in the p direction, E/N - Eq, only approach the ideal 
gas energy per particle in the p direction, hiVp, with de- 
creasing temperature as expected (see the arrow on the 
y-axis). Above Tc the total energies per particle for a 
= increase linearly with increasing temperature. How- 



ever, below Tc the total energy per particle decreases 
non-linearly with decreasing temperature, as expected 
theoretically Q. The specific heat can be calculated by 
differentiating the total energy per particle with respect 
to the temperature to locate the critical temperature Tc, 
where the specific heat has a peak in plotting the val- 
ues of the specific heat as a function of temperature. In 
the present calculation, however, it is not easy to see a 
peak due to the finite size of our system and large error 
bars. The scaled energies for a = lOaub lie well above 
the non-interacting curve both below and above Tc. The 
difference in the total energy per particle between for an 
ideal Bose gas and for a hard-sphere gas is visible at the 
low temperatures, and the clear difference increases as 
T decreases because of the short-range structure of the 
hard-sphere potential, i.e., the comparison between the 
total energy for an ideal Bose gas and that for a hard- 
sphere gas shows that the scaled energies per particle 
clearly shows the effect of the interaction strength. 

Finally, we calculated the superfluid fraction from 
PIMC data using the projected area |45|, not using the 
winding number because we used open boundary condi- 
tions along all three axies. The superfluid fraction de- 
pends on the rotation axis for the anisotropic trap. In 
particular, we calculated the superfluid fraction with re- 
spect to the symmetry axis, i.e., the z axis. Figure 4 
shows the superfluid fraction as a function of the aspect 
ratio A calculated using the PIMC method. In this cal- 
culation we used T/Tc — 0.4, N —27, and three different 
interaction strengths = 0, = anb, and as = lOajn, 
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quasi-lD gas shows no noticeable dependence. 

In order to compare our calculated PIMC data for the 
superfluid fraction with the theoretical value of an ideal 
gas |3^, we calculated the theoretical superfluid fraction 
for an ideal Bose gas using Eq. (14) of Ref. j3£]. The 
superfluid fraction Ns/N is 1 - Q/Qrig E3, where 6 is 
the quantum mechanical moment of inertia and Orig is 
the classical moment of inertia. In our case, we used 
the Tc with the first finite-size correction term (see, e.g., 
Eq. {m of Ref. 0). For this, we modified Eq^l4) of 
Ref. [33 by simply rewriting Eq. (14) of Ref. [33 using 
Eqs. (8) and (9) of Ref. Js^ to yield O 



FIG. 5: The calculated superfluid fractions along the z axis 
at the aspect ratio A =100 and N = 27 for three different 
interaction strengths = (circles), as — am, (diamonds), 
and fls = lOflfli, (triangles) as a function of the scaled tem- 
perature T /Tc- The inset shows the superfluid fraction along 
the z axis for a quasi-lD Bose gas (using the same symbols as 
in Fig. Ol. Solid lines show the theoretical superfluid fraction 
along the z axis for an ideal Bose gas calculated from Eq. ^ . 



(using the same symbols as in Fig. 3). For aspect ratios 
A < 100, the superfluid fraction decreases as the aspect 
ratio increases. However, the superfluid fractions do not 
depend on the aspect ratio within error bars for aspect 
ratios A > 100. Figure 4 also shows that the superfluid 
fraction decreases gradually as both the two-body inter- 
action strength and the aspect ratio increase. 

In Fig. 5 we present our calculated superfluid fractions 
for the aspect ratio A =100 as a function of the scaled 
temperature T/Tc- In the calculation, we used N= 27 
and three different interacting strengths as = 0, = aju,, 
and as = lOa^h (using the same symbols as in Fig. |2J|. 
The superfluid fraction decreases from unity at very low 
temperatures and becomes zero above the critical tem- 
perature Tc as the temperature increases. The superfluid 
fraction has a small value at the same temperature for 
an interacting Bose gas compared to a non-interacting 
Bose gas, which shows that the critical temperature for 
a repulsively interacting Bose gas is lowered compared to 
the ideal Bose gas as a result of interaction. The inset of 
Fig. 5 shows the superfluid fraction for a quasi-lD Bose 
gas 0]. 

The superfluid fraction for a quasi-lD Bose gas de- 
creases much more slowly than that for a quasi-2D Bose 
gas with increasing temperature. For example, at the 
scaled temperature, T/Tc = 1, the superfluid fraction 
for a quasi-lD system is significantly larger, about 0.65, 
compared to that for a quasi-2D gas, about 0.01. How- 
ever, the superfluid fraction for a quasi-2D gas shows a 
clear dependence on the interaction strength whereas a 



\N ) ^ 1 - (r/Tc)3 + imiiT/T^YikBTjnujpY ^ ^ 

In order to derive Eq. (14) of Ref. iSQi], Stringari used the 
semiclassical approximation in the so-called macroscopic 
limit. The solid lines in Fig. 5 and in the inset of Fig. 5 
show the resulting approximate superfluid fraction. For a 
quasi- ID Bose gas, the theoretical values agree well with 
PIMC data below T/Tc < 0.6 and after T/T^ = 0.6 the 
difference between the theoretical values and PIMC data 
gets larger as the scaled temperature increases. How- 
ever, for a quasi-2D system, the values calculated from 
Eq. JSJ lie well above our calculated PIMC data below Tc. 

In summary, using a finite-temperature path integral 
Monte Carlo technique we have investigated a trapped 
Bose gas in quasi-two dimensions interacting via a hard 
sphere potential whose core radius equals its correspond- 
ing scattering length. In order to enter a quasi-2D 
regime, we have changed the axial trapping frequency 
LUz ■ We have presented and analyzed our accurate PIMC 
results such as the total energy, the density profile, and 
the superfluid fraction as a function of temperature T/Tc 
at various values of the strength of the interaction and 
the aspect ratio A. We have compared our results with 
both the semi-classical approximation and the ideal Bose 
gas. We find that for repulsive interactions, the mini- 
mum value of the aspect ratio A, where the system start 
to behave quasi-two dimensionally, increases as the inter- 
action strength increases. In addition, for N — 27 and 
A = 100, the motion along the axial direction is indeed 
frozen out for all temperature range used in this paper 0.1 
< T/Tc < f -4 and the density profile n{z) becomes identi- 
cal to that of an ideal gas. Specifically, we have calculated 
the superfluid fraction using the projected area and have 
compared the superfluid fraction for a quasi-2D system 
to that for a quasi-lD Bose gas. The superfluid fraction 
for a quasi-2D Bose gas is distinctly different from that 
for both a quasi-lD Bose gas and a true 3D system, i.e., 
the superfluid fraction for a quasi-2D Bose gas decreases 
faster than that for a quasi-lD system and a true 3D sys- 
tem with increasing temperature, and shows a stronger 
dependence on the interaction strength. In addition, the 
superfluid fraction for a quasi-2D Bose gas lie well below 
the values calculated from the semi-classical approxima- 
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tion and the Kosterlitz-Thouless transition temperature 
decreases as the strength of the interaction increases. 
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